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Abstract. The main objective of this and its accompanying articles is to 
derive a mathematical theory associated with the thermohalinc circulations 
(THC). This article provides a general transition and stability theory for the 
Boussincsq system, governing the motion and states of the large-scale ocean 
circulation. First, it is shown that the first transition is either to multiple 
steady states or to oscillations (periodic solutions), determined by the sign of 
a nondimcnsional parameter K, depending on the geometry of the physical 
domain and the thermal and saline Raylcigh numbers. Second, for both the 
multiple equilibria and periodic solutions transitions, both Type-I (continuous) 
and Type-II (jump) transitions can occur, and precise criteria are derived in 
terms of two computable nondimcnsional parameters fei and 62- Associated 
with Type-II transitions are the hysteresis phenomena, and the physical reality 
is represented by either metastable states or by a local attractor away from the 
basic solution, showing more complex dynamical behavior. Third, a convection 
scale law is introduced, leading to an introduction of proper friction terms in 
the model in order to derive the correct circulation length scale. In particular, 
the dynamic transitions of the model with the derived friction terms suggest 
that the THC favors the continuous transitions to stable multiple equilibria. 
Applications of the theoretical analysis and results to different flow regimes 
will be explored in the accompanying articles. 



1. Introduction 

One of the primary goals in climate dynamics is to document, through careful 
theoretical and numerical studies, the presence of climate low frequency variability 
to verify the robustness of this variability's characteristics to changes in model pa- 
rameters, and to help explain its physical mechanisms. The thorough understanding 
of this variability is a challenging problem with important practical implications 
for geophysical efforts to quantify predictability, analyze error growth in dynamical 
models, and develop efficient forecast methods. 

Oceanic circulation is one of key sources of internal climate variability. One im- 
portant source of such variability is the thermohaline circulation (THC). Physically 
speaking, the buoyancy fluxes at the ocean surface give rise to gradients in temper- 
ature and salinity, which produce, in turn, density gradients. These gradients are, 
overall, sharper in the vertical than in the horizontal and are associated therefore 
with an overturning or THC. 
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The thcrmohalinc circulation is the global density-driven circulation of the oceans, 
which is so named because it involvse both heat, namely "thermo" , and salt, namely 
"haline". The two attributes, temperature and salinity, together determine the 
density of seawater, and the defferences in density between the water masses in the 
oceans cause the water to flow. 

The thermohaline circulation is also called the great ocean conveyer, the ocean 
conveyer belt, or the global conveyer belt. The great ocean conveyer produces the 
greatest oceanic current on the planet. It works in a fashion similar to a conveyer 
belt transporting enormous volume of cold, salty water from the North Atlantic to 
the North Pacific, and bringing warmer, fresher water in return. Figure [TTT] gives 
a simplified map of the great ocean conveyer and Figure |1.2| gives a diagram of 
oceanic currents of thermohaline circulation. 

In oceanography, the procedure of the Conveyer is usually described by starting 
with what happens in the North Atlantic, under and near the polar region see ice. 
There warm, salty water that has been northward transported from tropical regions 
is cooled to form frigid water in vast quatities, which results in a bigger density 
of seawater (unlike fresh water, saline water does not have a density maximum at 
4°C but gets denser as it cools all the way to its freezing point of approximatively 
— 1.8° C). When this seawater freezes, its salt is excluded (see ice contains almost 
no salt), increasing the salinity of the remaining, unfrozen water. This salinity 
makes the water denser again. The dense water then sinks into the deep basin of 
the sea to form the North Atlantic Deep Water (NADW), and it drives today's 
ocean thermohaline ciculation. 




Antarctica 



Figure 1.1. The diagram of the thermohaline circulation. The 
dotted line represents deep-water currents, while the solid line rep- 
resents shallow- water currents. 

The sinking, cold and salty water in the North Atlantic flows very slowly and 
southward into the deep abyssal plains of the Atlantic. The deep water moves then 
through the South Atlantic around South Africa where it is split into two routes: 
one into the Indian Ocean and one past Australia into the Pacific. As it contin- 
ues on its submarine migration, the current mixes with warmer fresh water, and 
slowly becomes warmer and fresher. Finally, in the North Pucific, the warmer and 
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deep - water currents 

Figure 1 . 2 . A schematic diagram of oceanic currents of the global 
thcrmohalinc circulation. 



fresher water upwells, while a shallow-water counter-current has been generated. 
This counter-current moves southward and westward, through the Indian Ocean, 
still heading west, and rounding southern Africa, then crosses through the South 
Atlantic, still on the surface (though it extends a kilometer and a half deep). It 
then moves up along the east coast of the North America, and on across to the 
coast of Scandinavia. When this warmer, less salty water reaches high northern 
latitudes, it chills, and naturally becomes North Atlantic Deep Water, completing 
its circuit. 

The THC varies on timescales of decades or longer, as far as we can tell from 
instrumental and paleoclimatic data [Martinson et al., 1995]. There have been 
extensive observational, physical and numerical studies. We refer the interested 
readers to [3J for an extensive review of the topics; see also among others l2Tfl 

HU EO EH EH Hfl 123 H3 111 U M E! ■ 

The main objective of this series of articles is to study dynamic stability and 
transitions in large scale ocean circulations associated with THC. A crucial start- 
ing point of this theory is that the complete set of transition states are described by 
a local attractor, rather than some steady states or periodic solutions or other type 
of orbits as part of this local attractor. Following this philosophy, the dynamic tran- 
sition theory is recently developed by the authors to identify the transition states 
and to classify them both dynamically and physically. The theory is motivated by 
phase transition problems in nonlinear sciences. Namely, the mathematical theory 
is developed under close links to the physics, and in return the theory is applied 
to the physical problems, although more applications are yet to be explored. With 
this theory, many long standing phase transition problems are either solved or be- 
come more accessible. In fact, the study of the underlying physical problems leads 
to a number of physical predictions. For example, the study of phase transitions 
of both liquid helium-3 and helium-4 [5TJ [SSJ HI] leads not only to a theoretical 
understanding of the phase transitions to superfluidity observed by experiments, 
but also to such physical predictions as the existence of a new superfiuid phase C 
for liquid helium-3. In return, these physical predictions provides new insights to 
both theoretical and experimental studies for the underlying physical problems. 
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The main objective of this article is to provide a theoretical framework for the 
general Boussinesq equations, which are basic equations describing the motion and 
states of large scale ocean circulations. Specific applications of the analysis and 
the results obtained in this article will be explored in the accompanying article in 
this series, including e.g. the effects of the turbulent frictions and earth's rotation, 
and circulations in different basins of the ocean. We remark that in collaboration 
with Hsia, the authors have done preliminary and related work on the dynamic 
bifurcations of the doubly-diffusive Boussinesq equations [121 [H] • 

To explain the main results obtained in this paper, we introduce a nondimcn- 
sional parameter K defined by: 



(1.1) K = sign(l - Le ) 



Le 2 ( I 



1 - Le V Pr 



1 + —)<7 C -R 



Here, as defined by (2.7 1, Pr is the Prandtl number, Le is the Lewis number, R is 



the thermal Raylcigh number, and R is the saline Raylcigh number. Moreover, the 



critical number a c is defined by (1.2 1, and can be casted as 

CT = min ^ 4 (j 2 Lr 2 + fc%- 2 + i) 3 = Afxr!+M^!±l)! 

c i 2L i 2+k2L * 2 jf^+kiL^ 

for some integer pair (ji,fci) such that j\ > 0, fei > 0, jf + k\ ^ 0. Here L\ 
and L2 are the nondimensional length scales in the zonal and meridional directions 
respectively (with the vertical length scaled to 1). 

The analysis in this article shows that the Boussinesq system undergoes a first 
transition either to multiple equilibria or to periodic solutions (oscillations), dic- 
tated by the sign of the nondimensional parameter K . In the case where K > O, 
the first dynamic transition of the system occurs as the R-Rayleigh number 

a = R-Le~ 1 R 

crosses the critical number a, leading to multiple equilibria. The transition can ei- 
ther be a Type-I (continuous) transition, or a Type-II (jump) transition, depending 
on the sign of the following parameter: 

1 — Le 2 

(1.3) bx = cr c ^^R. 

V ' Le 3 

The Type-II transition leads also to the existence of metastable stables, saddle-node 
bifurcations and the hysteresis associated with it. In this case, as the R-Rayleigh 
number goes beyond the critical value a c , the physical reality is represented by local 
attractors away from the basic equilibrium state. 

In the case where K < 0, the first transition of the Boussinesq system occurs as 
the critical C-Rayleigh number 

Pr +Le ~ 
77 = ~Pr~Tl~ 



crosses its first critical value 



(Pr+Lc )(1 + Le ) 



Vc = °c 

leading to periodic solutions (oscillatory mode). The spatiotemporal patterns de- 
picted by the oscillatory mode plays an important role in the study of climate 
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variabilities. As in the previous case, the transition can either by Type-I or Type- 
II, determined by the sign of another nondimensional parameter 62 defined by (4.2 1. 
As before, in the Type-II transition case, the physical transition states are more 
complex. 

The second part of this article addresses the discrepancies associated with the 
circulation scales. It is easy to observe that for large scale atmospheric and oceanic 
circulations, the critical convection scales are often to small. The effects of friction 
terms to the dynamics in geophysical flows have been studied by many authors; see 
among many others [2^1 121 [TOj |5TJ1 |5T] . 

The main objective of this part of the study is to propose a turbulent friction 
term in the Boussinesq equations, leading to correct convection scales. In particular, 
a basic convective scale law is derived in (5.20). With this scaling law, we are able 
to propose the correct friction terms in the Boussinesq equations to study THC. 
With this model, we show that the oceanic THC favors the continuous transition 
to stable multiple equilibria, rather than other types of transitions. 

This article is organized as follows. Section 2 introduces the classical Boussinesq 
equations and their nodimensional form. Linear analysis is given in Section 3, 
nonlinear dynamic transitions are studied in Section 4, concluding remarks are 
given in Section 5. The proof of the main theorems is provided in the appendix. 



2. Boussinesq equations 

As mentioned in the Introduction, to demonstrate the basic issues, we ignore 
the rotational effect of the earth, which will be studied in a forthcoming article. 
In addition, for simplicity, we regard the region of the Pacific, the Indian Ocean, 
and the Atlantic, where the ocean conveyer belt occupies, as a rectangular domain: 

= (0, ^i) x (0, 12) x (0, h), where l\ stands for the length of this conveyer belt, I2 
for the width, and h for the deep of the ocean. The motion and states of the large 
scale ocean are governed by the following Boussinesq equations (see among others 
[261 EQl Ell EH CU): 





du 


-(u. 


V)u = 


u/S.u (Vp + pgk) 

Po 


(2.1) 


or 

dS 

at" 


-(« 


V)T = 


k t AT, 




-(« 


V)S* = 


k s AS, 




divu 


= 0, 







where u — (1*1,1*2,143) is the velocity field, T is the temperature function, S is the 
salinity, k = (0, 0, 1), kt > is the thermal diffusivity, K5 > is the salt diffusivity, 
pa > is the fluid density at the lower surface X3 = for (x\, £2, £3) € f2, and p is 
the fluid density given by the equation of state 

(2.2) p = p [l-a(T-T ) + b(S-S Q )}. 

Here a and b are assumed to be positive constants. Moreover, the lower boundary 
(&3 = 0) is maintained at a constant temperature To and a constant salinity So, 
while the upper boundary = h) is maintained at a constant temperature Ti 
and a constant salinity Si- The case where T < T% and S > Si often leads to 
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thermal and solute stabilities, while each of the following three cases may give rises 
of instability: 



(2.3) 


T >T 1 , 


So < 5*i 


(2.4) 


r >r lj 


So > Si 


(2.5) 


T <T 1; 


So < Si 



The conditions (2.3 1- (2. 6) are satisfied respectively over some different regions of 



the oceans. In particular, in the high-latitude ocean regions the condition (2.3 1 or 



(2.4 1 is satisfied, and in the tropical ocean regions, the case (2.5 1 occurs. It is these 



properties ( 2.3 )-( 2.5 1 that give rise to the global thermohaline circulation, and the 



theoretic results derived in this section also support the viewpoint. 



The trivial steady state solution of ( 2. 1 )-( 2.2 1 is given by 



u u = 0, 

T° = T - (T - Ti)xa/h, 
S° = S - (S - Si)x 3 /h, 

P° = Po - 9Po[x 3 + ^(T ~ Ti)xj/h - |(Sb - Si)xl/h], 

where po is a constant. To make the equations nondimensional, we consider the 
perturbation of the soution from the trivial steady state 



u" = u- u°, 
S" = S- S°, 



rj-\ rpO 

P"=P-P°. 



Then we set 



x = hx' , 

u" — ktu' /h, 

S" = \So — S\\S' , 



t = h 2 t' /k t , 

T" = (Tq-T^T', 

p" = p Q vK T p'/h 2 . 



Omitting the primes, the equations (2.1 1 can be written as 
du 



at 

dT 

(2.6) Ik 
dS 



= Pr (Au - Vp) + Pr (RT - sign(5 - S l )RS)k - (u ■ V)u 
= AT + u 3 - (u ■ V)T, 
= Le A5 + sign(So - S^ua - (u ■ V)5 



at 

divu = 



for x — (x\, x%, X3) in the nondimensional domain 



fi= (O.Li) x (0,L 2 ) x (0,1), 
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where Li = li/h (i = 1,2), and the nondimcnsional parameters are given by 

R = — the thermal Rayleigh number, 

ktv 

R = — — ^ — the saline Rayleigh number, 

(2.7) k t v 
v 

Pr = — the Prandtl number, 

kt 

Le = — the Lewis number. 

Kj 1 

We consider the free boundary condition: 

du 2 du 3 OT 8S 

Ul = °> a — = a — = a — = a — = at Xl = 
ox i ox i ox i ox i 

,„„. „ dm du 3 dT dS n 

(2.8) u 2 = 0, -± = — ^ = — = — = at x 2 = 0, L 2 , 

OX2 OX2 OX2 OX2 

dui du-2 

u a = T = S = 0, — - = ~ = at x 3 = 0,1. 

The initial value conditions are given by 

(2.9) (u, T, S) = (u, f, S) at t = 0. 



For the problem ( |2.6[ )-(2.9 1, we set the spaces 

H={(u, T, S) e L 2 (ft) 5 I divu = 0, u • n\ an = 0}, 
Hi = {(u,T,S) G H 2 (n) 5 nH I T = S= at x 3 = 0,1 }. 
Let L\ = A + B\ : Hi ^ H and G . H\ ^ H be defined by 
Aip = P(Pr Au, AT, Le AS) , 

= P(Pr (RT - i?5sign(5 - Si))fc, u 3 , sign(5 - £1)1*3) 
G(^) = P((«-V)«,(t*-V)T, («-V)5), 

for V = (u,T,S) € Hx,\= (R,R). Here P : L 2 (n,R 5 ) -> iJ is the Leray projec- 
tion. Then the problem ( 2.6 1-( 2.9 ) can be written as 

(2.10) d i=L^ + G W , 

3. Linear Theory 

3.1. Eigen-Analysis. To understand the dynamic transitions of the problem, we 
need to study the following eigenvalue problem for the linearized equations of (2.6)- 
@: 

Pr (Au - Vp) + Pr (RT - RSsign(S - Si))fc = 0u, 

(3.1) AT + . 3 =/?T, 

Le AS + sign(5 - Si)u 3 = f3S, 

divu = 0, 



supplemented with (2.. 
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We proceed with the separation of variables. By the boundary condition (2.8 1, 
ip — (it, T, S) can be expressed in the following form 

ui — Ujk(xs) sin jaiirxi cos kot2'KX2, 

U2 = Vjk{xz) COS jOliTTXi sin kcC2lTX2 , 

M3 = Wjk{x^) cos jaiirxi cos kot2iTX2, 
T = Tjk{xz) cos jaiirxi cos ka2irX2, 
S = Sjk(xs) cos jaiirxi cos ka2TTX2, 
P = Pjk{x^) cos jaiirxi coska2^X2, 



for integers j and fe, where ccj = L i (« = 1, 2). Plugging (3.2 1 into (3.1 1, we obtain 
the following sets of ordinary differential equation systems: 

Pr D jk u jk + Pr jaxirpjk = (3u jk , 

Pr D jk v jk +Pr ka 2 Trpjk = 0Vj k , 

Pr D Jfc u; ifc - Pr DP jk + Pr i?T jfc - Pr Rsign(S - S 1 )S jk = (3W jk , 

(3.3) D jk T + w jk = /3Tjk, 

LeDjkSjk + sign(S'o - Si)w jk = PS jk 

jaiTTU jk + ka 2 TTv jk + Dw ]k = 0, 

Dujk = Dvjh = Wjk = T jk = S 3 k =0 at x 3 = 0, 1, 

where 

d d 2 

7-1 a 7-1 u 2 2 2/ -2 2 i 7 2 2\ 

£>=— , D ]k = — -a ]k , a jk = TT [j a x + k a 2 ). 
If Wjfc 7^ 0, equations (3.3 1 can be reduced to a single equation 
[(D jk - /3)(Le D jk - /3)(Pr £ jfc - /9)D ifc 

(3.4) + Pr i?a 2 fe (Le D jk - 0) - Pr Ra 2 k (D jk - (3)}w jh = 0, 
Wjfe = D 2 w jk = D 4 w jk = D 6 w jk = at x 3 — 0, 1, 



It is clear that the solutions of (3.4) are sine functions 
(3.5) w jk = smlirx 3 for (j, k, I) £ I? x N. 



Substituting (3.5 1 into (3.4 1, we sec that the corresponding eigenvalues (3 of Problem 
(3.11 satisfy the cubic equation 

(3.6) 7 ] kl (j^ + P) (Le 7 ? w + (3) (Pr j 2 kl + 0) 

-Pr Ra%(Le ~? 2 kl + 0) + Pr Ra 2 jk (~/jki + 0)=O V(j, fc, Z) G Z 2 x N. 

where 7 2 fc; = a| fc + Z 2 tt 2 . 

Moreover, to determine Uj k {x 3 ),Vj k {x 3 ),Tj k {x 3 ), and Sj k {x 3 ), we derive from 



(|3.3] ) that 

(3.7) (Pr £> jfc - ,3)u jfc = -^(Pr £> jfc - /8)Du; Jfc> 

(3.8) (Pr D jk - 0)v jk = -^^(Pr D ifc - /SJDwj*, 

(3.9) {D jk - 0)T jk = -w Jk , 

(3.10) (Le Djk - (3)Sjk = -sign(5 - SJwjk 
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With the above calculations, all eigenvalues and eigenvectors (3.1 1 can be derived 
and are given in the following three groups: 

1. For (j, k, I) — (J, k, 0), we have 

Pjko = -Pr a% = -Pr {fa 2 + k 2 a 2 2 ), 

ipjko = (ka 2 sinjaiTTXi cos ka 2 irx2, —ja± cos jaiirxi sin kai~KXi^ 0, 0, 0) 

2. For (j, k, I) = (0, 0, I) with I ^ 0, we have 

Pool = -' 2 7r 2 , Pool = ~ Le ' 2?r2 ' 

"0oo; = (0,0,0, sin Z7ra 3 ,0), ^ooz = (0, 0, 0, 0, sin^Tra^). 



For general (j, k, I) with j 2 + k 2 ^ and I ^ 0, the solutio ns /? of (3.6 1 are 
eigenvalues of (3.1 1. Let /3j- fe; (1 < i < 3) be three zeros of (3.6 1, with 

Re[3} kl >Re(3 2 kl >Refi kl . 



Then, by (3.2 1,(3.5 1, and ( 3.7 1-( 3.10 1, the eigenvector ip* kl corresponding 
to (3j kl can be written as 



(3.11) 



jkl 



= < 



Uj k i sin jaiirxi cos ka2tTX2 cos liTX3 
Vj k i cos jaiirxi sin ka2TX2 cos Ittx^ 
Wj kl cos jaiirxi cos fca27ra:2 sin/Va^ 
Tj fe ; cos jaiirxi cos ka 2 irX2 sin ?7ra;3 
_ S™ fe ; cos jaiirxi cos £;a27ra;2 sin /7ra3, 



where 



= 1, 



T. 



jctilir 2 
1 



jail 



j 2 a'f + k 2 a 2 



V jkl = 



ka 2 lTT z 

^2 



ka2l 



a 



jk 



j 2 af + k 2 a\ 



jki 



'y 2 + B 2 ' 

l jkl ' > J jkl 



Vjkl 



signup - Si) 
^1%1 + P 2 



jkl 



The adjoint equations of (3.1 1 are given by 



(3.12) 



Pr (Au* -Vp*) + (T* + sign^o - Si)S*)k = 0u* 

AT* +Pr Ru* 3 = j3T*, 

Le AS* - Pr 7?sign(5 - Si)u* 3 = (3T* 

divw* = 0. 



Thus, the conjugate eigenvector ip^* kl of (3.121 corresponding to fi^ kl are as follows 

i l * k i sinjaiTrxi cos ka 2 irx2 cos Ittx^ 



(3.13) 



jkl 



jkl 

A* 
"jkl 



cos jaiirxi sin ka2^x 2 cos lirx 3 
cos jaiirxi sin ka.i'KXi sin (Va^ 
TJy cos jaynxi sin ka2^x 2 sin ^£3 
t Sj^j cos jaiirxi sin fca27r^3 sin ^23, 
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where 



l jkl 



u 3 kl - - 



y 2 



u jki 



fca2^7r 2 



a' 



jki 



T 



Pr R 



jkl 



S 



J jki 



jki 



-sign(5 - 5i)Pr R 



Rej 2 



jki 



■Pi 



k-i 



Thus, all eigenvectors of (3.1) consist of ipjko, V'ooi' V'ooi' an( ^ V^/d 



conjugate eigenvectors of (3.12) consist of tp* k0 
and il)f kl as in (Ell. 



(i = 1,2). All 

002' 



3.2. Principle of exchange of stabilities (PES). The linear stability and in- 
stability are precisely determined by the the critical-crossing of the first eigenvalues 
of (3.1l, which is often called PES. For this purpose, we only need to study the 
solutions (3 of (3.6), which is equivalent to the following form 

(3.14) 



jkl 



/3 3 + (Pr + Lc + l) 7 | fei /3 2 + [(Pr + Le + Pr Le )~f% 
- Pl ' <x%7iki(R ~ R)]P + Pr Le j% - Pr a| fe (Le R - R) = 0. 



As we shall see, both real and complex eigenvalues can occur. The real eigenval- 
ues often lead to transition to steady state solutions, and the complex eigenvalues 
gives rise oscillations. As we mentioned in the INtroduction, these solutions are 
related to low frequency variabilities of the oceanic system. 

First, we discuss the critical-crossing of the real eigenvalues. To this end we need 
to introduce a nondimensional parameter, called the R-Rayleigh number, defined 
by 

gh 3 



(3.15) o = R- Le _1 # = — (a(T - T t ) - 6Lc _1 (5 

ktv 



Si))- 



It is clear that (3 — is a zero of ( |3.14 1 if and only if 
Le 1 % l -a%(LeR-R) = 0. 



In this case, we have 
(3.16) 



-v 6 

ljkl 



> G r 



a 



jk 



where a c is defined by (1.2 1, and is called the critical R-Rayleigh number. 

Next, we consider the critical-crossing of the complex eigenvalues. For this case, 
we introduce another nondimensional parameter, called the C-Rayleigh number, 
defined by 



(3.17) 



r) = R- 



Pr +Le 
Pr + 1 



R 



a(T Q - Ti) 



(Pr +Le)6 
Pr + 1 



(So - Si) 



Let ipo (po 7^ 0) be a zero of (3.141. Then have 

pi = (Pr + Lc + Pr Lc ) 7 | fci - Pr a 2 jkl j*{R - R), 



Pl = 



Pr Lc 7 L-Pr a 2 (LeR-R) 



'jkl 

7 2 fei (Pr +Lc 



1) 
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Therefore, equation (3.141 has a pair of purely imaginary solutions ±«po if and only 

Pr a%(Pr + Le + 1)(R-R) 



Pr Le h% 



if the following condition holds true 
(3.18) (Pr +Le + l)(Pr + Le 

= Pr Le - Pr a%(Lc R - R) > 0. 
It follows from (3.18) that 

R 



Pr 



V 



Lg -r 



(Pr +Le ){l + Le)Hju 



Pr + 1 Pr 
Hence we define the critical C-Rayleigh number by 



a 



(3.19) 



mm 

U,k t i)ei 



(Pr + Le )(1 + Le ) l% t _ (Pr + Le )(1 + Le ) 7, lfcl i 



Pr 



Pr 



for the same integer pair (jx,hi) as in (1.2 1. 



Definition 3.1. Let (ji,kx) satisfy {1.2). We call a c (resp. rj c ) the first critical 
Rayleigh number if for all eigenvalues [3 of \3. i| ) we have 

Rej3(a c ) < (resp. Ref3{r) c ) < 0). 

The following theorem provides a criterion to determine the first critical Rayleigh 
number for o~ c and rj Cl and the proof is given in the appendix. 

~,2 



Theorem 3.1. Lethe ^ 1, and (ji,ki) satisfy (1.2), and 7 = 7^.^, a 
Then we have the following assertions: 



(1) If K > where K is defined by (1.1), then the number o~ c is the first critical 
Rayleigh number, and 



(3-20) p) ikil 



< if a < a Cl 
= ifa = a c , 
>0 ifa>a c , 



(3.21) Refi M (* c ) < 



for all (j,k,l) not satisfying (1.2 1. 



(2) If K < 0, then the number r\ c is the first critical Rayleigh number, and 

<0 ifv<Vc 

= ifr) = T) c , 

>0 ifn>j] c , 



(3.22) Refi) ikil = ite$ ifcil 



for (j,k,l) not satisfying (1.2) 



(3.23) Ref3V kl (r) c ) < 

This theorem provides a precise criteria on if the first unstable mode corresponds 
either to multiple equilibrium modes or to oscillatory modes, corresponding to 
steady state patterns or to spatiotemporal patterns. 

4. Nonlinear Dynamic Transitions 
4.1. Transitions to multiple equilibria. By Theorem |3.1[ we know that under 



the conditions K > with K defined by (1.1 1, the Boussinesq equations (2.6 1 with 
(2.8 1 will have a transition at a = o~ c from real eigenvalues. In this section, we 
study the transition from a real simple eigenvalue. We know that generically, the 
first eigenvalues of (3.1) are simple. Hence we always assume that the first real 
eigenvalue of (3.1 1 near a = a c is simple. Then we have the following results. 
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Theorem 4.1. At 



that Le ^ 1, K > 0, and fei > 0, where K is defined 



by (1.1), and b\ is defined by (1.3). Then the problem (2.6) with (2.8) under- 



goes a Type-I (continuous) transition at the critical R-Rayleigh number a c , and the 
following assertions hold true: 

(1) If the R-Rayleigh number a crosses a c , the problem bifurcates to two steady 



state solutions t/>f = «, Tf , Sf ), i = 1,2. 

(2) There is an open set U C H with ip = G E7 snc/i £/ia£ {/ 
I7i n E/ 2 = 0, V> = G dUi n <9f7 2 , awd ^z " G ^ attracts [/< (i = 

(3) i/jf can oe expressed as 

C = (-l) l av ^)^ lfcll +o(/3 1/2 (a)) 1,2, 



= f/i 
1,2). 



U 2 



where a > is a constant, (3{u) = /3j lfcll (c) is the first eigenvalue satisfying 



(4) 



(B.2), and iph^i is the first eigenvector given by (3.11). 
If ipQ G Ui (i = l,2j ; i/iere exists a to > such that when t > to, the 
velocity component ofipit, tpo) is topologically equivalent to the structure as 
shown in Figure \6~1\ with either the same or the reversed orientation. Here 
ip(t,ipo) is the solution of \2.6\j with (2.8) with initial data ipo- 



Theorem 4.2. Assume thatLe ^ 1, K > 0, and b\ < 0. Then the problem (2.6)- 
(2.8) undergoes a Type-II (jump) transition at a c , and the following assertions hold 
true: 

(1) The transition of this problem at a c is a subcritical bifurcation, i.e. there 
are steady state solutions bifurcated on a < a c , which are repellors, and no 
steady state solutions bifurcated on a > a c . 

(2) There is a saddle-node bifurcation of steady state solutions from (ipl,a*) 
and with a* < o~ c , and there are four branches of steady state 
solutions < * < 4) in which and tp^ are bifurcated from (ip*,a*), 
and tp3,ip0 bifurcated from (ip2,o~*) on a > a*, such that 

lim iPl= lim V'=0, <,<^0, 

a — >a c a — >a c 

and and ^ are attractors for a* < a < a c + e with some e > 0. 
4.2. Transition to oscillatory spatiotemporal patterns. If the condition K < 



holds true, by Theorem 3.1 the problem (2.6 (-(2.8 1 will have a transition to 



the periodic solutions at the critical Rayleigh number r\ c given by (3.191. In this 



section, we discuss the transition from complex eigenvalues. It is easy to see that 
as Li/L 2 ^> 1,77c is given by 



Vc 



(Pr +Lc )(1+Le )7j l0 i 



Pr 



a h0 



(4.1) 

for some (j'i , ki) 
hold true. 

In the following, we introduce a parameter which provides a criterion to deter- 
mine the transition type: 



(jijO) with ji > 0. Here, we always assume the condition (4.1) 



(4.2) 



R c (R2Bi — 12-62)63 + i? c (/2-E>i + i? 2 B 2 )C , 4 
+i?(i? 3 5i - I 3 B 2 )C 5 + R(I 3 Bt + R 3 B 2 )C e , 



DYNAMIC TRANSITION THEORY FOR THERMOHALINE CIRCULATION 



13 



where a 2 = a 2 i0 = j 2 a 2 Tr 2 , and 

R 2 ^ReA 2 (p 1 j2 (R C2 )) 

R 3 = ReA 3 



(a 2 + it 2 ) 2 + p 2 ' 
Le (a 2 + 7T 2 ) 



Le (a 2 + 7T 2 ) 2 + p 2 
I 2 =ImA 2 (pUR C2 )) = 



I 3 = ImA 3 = 



(a 2 + it 2 ) 2 + p 2 ' 
P 



Le 2 (a 2 + 7r 2 ) 2 + p 2 ' 
B\ = RI 3 — R C I 2 , 
B 2 = RI3R3 — RI 2 R 2 , 

3i? 2 , P 2 R2 P^h 



C 3 
C 4 

c 5 
c 6 



2tt 27r(47r 4 + p 2 ) 4tt 4 + p 2 

2 



h P 2 h pitR 



2ir 27r(47r 4 + p 2 ) 4tt 4 + p 2 ' 
3i? 3 P 2 R 3 ph 



27rLe 27rLe (4?r 4 Le 2 + p 2 ) 47r 4 Le 2 + p 2 ' 

h P 2 h pRs 

2nLe 27rLe(47T 4 Le 2 + p 2 ) 47r 4 Le 2 + p 2 ' 



2 Pr(l-Lc)fl a 2 2 2 22 

Pr + 1 (ar + tt z ) 
(Pr + Le )(1 + Le ) (a 2 + tt 2 ) 3 Pr + Le - 
Pr a 2 + Pr + 1 

The following theorem characterizes the transition to periodic solutions and the 
type of the transition. 



Theorem 4.3. Let b 2 be given by (4-2), and K < with K defined by (1.1). Then 
for problem (2.6)-(2.8), we have the following assertions: 

(1) Ifb 2 < 0, then the problem has a Type-I transition at n — n c . and bifurcates 
to a periodic solution on n > n c which is an attractor, the periodic solution 
can be expressed as 

ip = x(t)Reip jl01 + y(t)lmtp jl01 + o(\x\ + \y\), 

()= l N ) p ' 

where I3{rf) = \(r)) + ip(rj) is the first eigenvalue with X(rj c ) = Re(3(j] c ) = 
0,p(r] c ) — p as given in (4-2). 

(2) Ifb 2 > 0, then the transition at r\ = r\ c is of Type-II, and there is a singular 
separation of periodic solutions at some (T*,rf) for rf < i] c , where T* is a 
periodic solution. In particular, there is a branch T n of periodic solutions 
separated from T* which are repellors, such that r„ — > as ?; — > n c . 
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The following example shows that both cases with b 2 < and 62 > may appear 
in different physical regimes. 

Let p 2 = 0. In this case we have 

(Pr + l)Le 2 7 6 f > 



R 



Rr = 



Pr (1 - Le ) a 2 { < 

Pr +Lc 7 6 
Pr (1 - Le ) a 2 ' 
Thus, the number 62 defined by (4.2) is given by 

b 2 = - 7 - 12 (i? c -Lc" 4 fl)(i? c -Le- 3 J R) 
-1 

Pra 2 (l-Le)Le 
It is clear that 



for 1 > Le 
for 1 < Le 



3 (Le (Pr + Le ) - Pr - l)(Le 2 (Pr + Le ) - Pr - 1). 



< 
> 



if Le < 1, 
if Le > 1. 



5. Convections Scales and Dynamic Transition 

5.1. Convection scale theory. We investigate a phenomenon appearing in con- 
vection problems. For this purpose, we start with the classical Rayleigh-Benard 
convection as studied in |16[ 119] . In their nondimensional form, the Boussinesq 
equations are given by 



du 



dt 



(5.1) 



1 

Pr" 

dT 

~dt 
divu = 0. 



+ (u ■ V)u + Vp 



Au - \fRTk = 0. 



(u ■ V)T - \fRu3 — AT = 0, 



With the free boundary condition 



(5.2) 



u n = 0, 

T = 
dT 



On 



= on dfl, 



On 



= 



at X3 = 0, 1, 

at x\ — 0, L\ or x 2 = 0, L 2 . 



the critical Rayleigh number is 
(5.3) 



P -min (Q2 + 7r2)3 - (Qg + 7r2)3 - 277f4 
iXr — mm — — , 



where a 2 ~ and for the Dirichlet (rigid) boundary condition, by Chandrasekha 
H, R c is 

(5.4) R c = 1700. 



By (5.3 1 and (5.4 1, as the horizontal scale L of the fluid is much larger than its 

convection appears at the tempe: 

for the free boundary condition, 



height h, i.e. L 3> h, the thermal convection appears at the temperature difference 

27KI/7T 4 

_ J Agah 3 



(5.5) 



T0-T1 



gah 3 



for the Dirichlet boundary condition. 
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From (5.5) we see that the critical temperature difference satisfies 

(5.6) AT C = T - Ti -> oo if h 

(5.7) AT C = T Q - Ti -> if ft 



0. 



oo. 



Physically, (5.6 1 amounts to saying that for a given fluid there is a minimal size 
/lo > such that when h < ho, no convection occurs for any AT = To — T\ > 0. 
This is in agreement with the fact that in the real world, when h is small, it is hard 
to maintain a high temperature gradient to generate the vertical convection. 

However, (5.7 1 has certain physical discrepancies with the energy conservation 
of the system. In fact, as h increases, it needs more energy to overcome the friction 
force to drive fluid to move. With AT representing this driving force, (5.7) shows 
that when h increases, AT decreases. 

Here we present a method to resolve this issue. We note that the critical value 
a 2 = 7r 2 / L 2 stands for the horizontal convection scale L (h as its unit): 

L = 7r I a . 



Hence, at the critical value a 2 
(5.8) 



\ in (5.3 1, the convective cell size is 



When h is not large, this value (5.8 1 is a good approximation, we need to modify 



the model to accommodate the case when h is large. 

Instead of (5.8), we propose the convection scale L c to be an increasing function 
of h: 



(5.9) 



Li 



ip(h), such that -0(0) = 2, ip'(h) > 0. 



Thus, the critical Raylcigh number R c reaches its minimal at 



(5.10) 

rather than at a 2 = it 2 /2. 

The convection scale law of (5.9 1 and (5.101 should be reflected in the hydro- 
dynamical equations. In other words, we need to revise the Boussinesq equations 
leading to new critical values given by (5.9 1 and (5.10). The method we propose 
here is to include a (turbulent) friction term: 

(5.11) F = -(<7 Ul,<JoU 2 ,<7lU 3 ) 

in the dimensional form of the Boussinesq equations, where <7j depends only on h 
satisfying 

(5.12) cr,(0) = 0, o--(/i)>0, i = 0,l. 



Physically, (5.111 and (5.121 stand for the added resistance generated by pressure, 
also called the damping terms. One can also argue this is due to the averaging to 
resolve small scale eddies. 

In their nondimensional form, the revised Boussinesq equations read 

du 
~di 

(5.13) dT 
dt 

divu = 0, 



+ (u ■ V)u = Pr [Au 
+ (u ■ V)T = AT 



Vp-f + RTk], 



"3, 
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where / = (5qUi, 5qU2, 81U3), and 



(5.14) 



Si{h) = — Oi{h\ (i = 0, 1). 



For the revised equations (5.131 with the free boundary conditions (5.2), the 
critical Rayleigh number is 



(5.15) 



Rc 



(a 2 + ir 2 )5 1 + 



(a 2 + ir 2 ) 3 + 7r 2 (a 2 + 7r 2 )5 



Let g(x) = (x + tt 2 )Si + ((x + n 2 ) 3 + tt 2 (x + tt 2 )8q)/x. Then, x c — a 2 satisfies 
g'(x) — 0. Thus we get 



(5.16) 



Sxal - (ai + tx z Y - TT z (a 2 c + tt z )S + 3a z c (ai + tt 2 ) 2 + ir z a% = 0. 



We assume that 81 > S for h + 0. By a 2 = n 2 /L 2 from ( |5.1l[ ), ( |5.12| ) and ( |5.16[ ) 
we infer the convection scale law of (5.9) and (5.101. In particular, when h is large, 
<5i ^> 1 and 5i 3> ^o, from (5.151 and (5.161 we derive that 

(5.17) 



(5.18) 
(5.19) 



4 _ TT^^+So) 

a ^ — * — ' 

L 2 c =5y 2 /(n 2 + S )^ 2 , 
R c ^ tt 2 <Si + ir 2 5 L 2 c , 



These formulas (5. 17 1- (5.191 are very useful in studying large-scale convection mo- 
tion for the height of fluid h > 10m. It can fairly solve the contradiction caused by 
(5.7 1. In fact, by (5.141, (5.18) and ( 5.19[ ) we get the critical temperature defference 

ir 2 Kai(h) 



gah 



for h ^> lm. 



It suggests that ai(h) = c ■ h 2 , leading to the critical temperature gradient being 
independent of vertical scale h. 

In summary, by (5.141, we propose to take 

(5.20) 5 = Coh^/v, 61 = dh^/v, 
where Co and C\ are constants. 

5.2. Revised Boussinesq equations of the large scale ocean and critical 
parameters. Consider the equations (2.1l with a damping term 5u. In their nondi- 
mensional form, the revised equations of (2.1 ) are given by 

DT 

(5.21) 9i 

as 
at 

divu = 



Au-5u + RTk - sign(5 - S^RSk -Vp - (u • V)u, 
= AT + u 3 - (u.V)T, 
= Le AS + sign(S - Si)u 3 - (« ■ V)S, 



supplemented with the boundary condition (2.8), where Su — (SqUi, SqU2, S1U3). 
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As in (3.6), the eigenvalues of the revised problem satisfy the following equation 
(5.22) a%(^k + /3)(Le j] h + 0)(Pr 7 2 fe + Pr S 1 + (3) 



+n 2 (j 2 k + (3)(Le 7 2 fe + /3)(Pr 7 2 , + Pr S + /?) 
Pr Ra%(Le 7 2 fe + (3) + Pr Ra%( 7 % +0)=O, 



which is equivalent to 
(5.23) 

7 2 ^ 3 + [ 7 4 (1 + Pr + Lc ) + Pr (a 2 ^ + 7r 2 <5 )]/? 2 

+ [(Pr + Lc + Pr Le ) 7 6 + Pr (1 + Lc ) 7 2 (a 2 <Si + tt 2 5 ) - Pr a 2 (R - R)}[3 
+ 7 8 Pr Le + 7 4 Pr Le {a 2 5 x + 7r 2 <5 ) - Pr 7 2 a 2 (Le R-R) = 0. 



Let (3 — 0, then we derive from (5.23 1 the critical Raylcigh number for the first real 
eigenvalues as 



(5.24) 



Rc 



Le R + min 



7 fli + 2 



where 7 2 = a 2 + 7r 2 . Let /3 = ipo (po ^ 0) in (5.23), then we get the critical 
Rayleigh number for the first complex eigenvalues as follows 

(5.25) r C2 =min(r 1 + r 2 + r 3 ) ! 

a 2 

where 

= (Pr + Lc ) 7 4 + Pr (a 2 S 1 + rr 2 ^) - 
1 (Pr +1) 7 4 + Pr {a 2 5 1 + tt 2 6 ) ' 

(Pr + l)(Pr + Lc )(1 + Lc ) 7 10 + Pr (1 + Lc )(1 + 2Pr + Lc )^{a 2 5 l 



7T 2 «5 ) 



Pr a 2 ((Pr + 1) 7 4 + Pr {a 2 5 1 + tt 2 S )) 
Pr 2 (l + Le )7 2 (a 2 ^ + 7r 2 <5 ) 2 



' Pr a 2 ((Pr + 1) 7 4 + Pr {a 2 5 1 + 7r 2 <5 )) ' 

It is readily to check that the first eigenvectors of the revised problem are the 



same as that in (3.11) and (3.131 



For (5.211, we need to modify the transition theorems, Theorems 4.1 4.2 and 



4.3 i.e., to modify the numbers &i and b 2 given by (1.3) and (4.2 1 associated with 



the new parameters (5.241 and (5.25). Thus, we can make a comparision between 



the new and old theories. 

5.3. Revised dynamic transition theory. Linking to the large scale oceanic 
circulation, we now study the dynamic transition problem for the revised equations 



(5.21 1. We start with the introduction of a few typical physical parameters. 



Physical Parameters in Oceanography. In oceanography, main physical pa- 
rameters are listed as 



(5.26) 



Pr 



Le = 1(T 



1.1 x 10~ 6 m 2 s-\ 



1.4 x 10~ 7 m 2 s~\ 



a 
h 



2.1 x 10" 4 K~ 
: 4 x 10 3 m. 



The haline contaction coefficient b satisfies 

^I^xKT 3 , 
p as 

P= Sp s + (1 - S)p w , 
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where S is the concentration of salt with unit 0/qq(j>su), ps and pw are the densities 
of salt and water given by 



p w = 10 3 kg/m 3 , p s = 1.95 x 10 3 kgmT 3 . 

At the sea water density p = 1.01 ~ 1.05 x 10 3 kgmT 3 , the halinc concentration 
coefficient b is about 

b = \p s - p w ) x 10~ 3 0.92 x lQ- 3 (vsu)- 1 . 
P 

Thus, for the oceanic motion we obtain the thermal and saline Rayleigh numbers 
as follows 



(5.27) 
(5.28) 



R = ga{T ° Tl) h 3 = 0.86 x 10 21 (T - T^C- 1 ], 



R 



gb{So - Si) 3 

KV 



h 3 = 3.75 x 10 21 (S - S^ipsu)- 1 . 



For the thermohaline circulation, we know that its scale is tens of thouthands 
kilometer, namely 



(5.29) 



L c = 0(10 4 ) (h = 4km as unit). 



For the damping coefficients Sq and Si, by (5.20), we have 



Si = C l h i /u = 2.33 x 10 20 Ci m 2 



= 0,1). 



Phcnomcnologically, the constants C, depend on the density p, or equivalently on 
the pressure p, and we assume that C\ is proportional to p. Note that the ratio of 
densities of water and air is about 10 3 , i.e. 

Pw/Pa = 10 3 . 

Therefore, the constant Ci of water is about 10 3 times that of air. Thus, by (5.201, 
for water we take 

(5.30) 



Ci w 10 3 (m 2 • s)- 1 . 



However, the constant Co depends also on the smoothness of the lower surface. 
Because the bottom of the sea is more smooth than the lands, this constant Co of 
sea water is not larger than that of air. Thus, we take 



(5.31) 



Co 



1 



x 10 



-12 



Then, we obtain the damping coefficients as 

(5.32) 5 = 1.17 x 10 s , Si = 2.33 x 10 23 . 



We remark that (5.301 and (5.311 are theoretically estimated values. However, 



they do provide an fairly reasonable explanation to the critical Rayleigh numbers 
and the convective scales in both the atmospheric and oceanic circulations. 



Critical Rayleigh numbers. By (5.261 and (5.321, Ti ~ R, where Ti is as 



defined in (5.251. Thus, the critical Rayleigh numbers (5.241 and (5.251 can be 
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approximativcly expressed as 
(5.33) R C1 - Lc^R = n 



(tt 2 + q 2 )S 1 + 



7r 2 (7r 2 + a 2 )S 



(5.34) R C2 - AS (1 + Le ) min 

a 2 

= ( 1 + Le ) min 

a 2 

Thus, they have the same a 2 given by 
(5.35) 

The convective scale is 



Pr+1 (a 2 + 7T 2 ) 5 2Pr+l (a 2 + tt 2 ) 3 



Pr a 2 (Si + tt 2 S Pr 



+(a 2 + tt 2 )5 1 + 
n 4 S 



n 2 (a 2 + n 2 )S 



(a 2 + tt 2 )Si + 



2 2 
a„ = 7T 



1/2 



2.24 x 1CT 7 



(5.36) 



1/4 



0.67 x 10 4 . 



The critical Rayleigh numbers for real and complex eigenvalues are respectively 
given by 



a c = (n 2 + a^S, + ^ = 2.33 x 10 23 , 

at 



(5.37) 

(5.38) 77 C = (1 + Lc )a c . 

Here, as in (3.15) and ( 3.16[ ), by (5.331 and (5.34), we define the R-Rayleigh number 
a and C- Rayleigh number rj as follows 



(5.39) 



<r = R-Le~ 1 R, i] = R-R. 



It is seen that the theoretical value (5.36) agrees with the realistic length scale 
given in (5.29 1. 



Revised transition results. From ( 5.37 1- ( 5.39 1 we can see that if 

a = R — Lc _1 i? = (T c , and rj = R — R < rj c — (1 + Le )er c , 

then R < Le 2 a c /(1 + Le ). In this case a c is the first critical Rayleigh number. 
In addition, if 

R-Lc~ 1 R= (1+Le )a c , R-R><r c , 

then R > Le 2 cr c /(1 + Le ), and r\ c is the first critical Rayleigh number. Thus, by 
Theorem 3.1 we obtain the following physical conclusions: 

Physical Conclusion 5.1. For the thermohaline circulation, we have the following 
assertions: 



(1) // 



R < Lc z <7 c = 2.35 x 10 19 , 



the number a c given by (5.31) is the first critical Rayleigh number. 
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(2) // 



R > 



Le 2 



1 -Lc 



-o> = 2.35 x 10 



1!) 



the number r\ c given by (5.38) is the first critical Rayleigh number. 



Note that the first eigenvectors of the linearized equations of (5.21 1 are the same 



b\ should be as in (1.3 1; namely 



as that of (3.1|. For the transition of (5.211 from real eigenvalues, the parameter 



h = <t c - 



Le J 



Here R — a c 4 
Thus we have 

(5.40) 



Le _1 /i, Le 



io- 



cr c = 2.33 x 10 2d , Pr 



b x 2.33 x 10 1 



where R is give n by (|5.28|) 
By Theorem 



5.3 



and of. = 2.24 x 10" 



R, 



as R < 2.35 x 10 19 the transition of (5.211 is at a = a c . The 



following theorem is a revised version of Theorems 4.1 and |4.2 



Physical Conclusion 5.2. Let R < 2.35 x 10 19 , and b\ be as in (5.40). Then 
the problem (5.21) with (2.8) undergoes a dynamic transition from a — cr c , and the 



following assertions hold true: 

(1) If R< 2.33 x 10 17 , the transition is continuous, and the problem bifurcates 



to two steady state solutions ip* and ip% f or a > °~c, a as in (5.39), ipl and 
ip2 have the expressions as in Assertion (3) of Theorem J^.l, with U\ and 
U% as their basin of attractions in H respectively. 

(2) 7/2.33 x 10 17 < R < 2.35 x 10 19 , the transition is jump, and there are two 
saddle-node bifurcations from (ipij°~*) an d (V'Si' 7 *) with a* < a c . 

(3) When the transition is Type-I, if the initial value tj) € Ui, then there is 
a time t such that as t > t , the velocity component u, in the solution 
ipit^ipo) with initial value ipg, is topologically equivalent to the structure as 
shown in Figure \6.1\ 



For the transition to periodic solutions, the revised parameters are as follows 
a 2 = a\ = 2.24 x 10~ 7 < 1, 

2 2 i 2^2 

7 = a + 7T = 7T , 

R = R C2 =R+T) C (r] c as in (10.208)). 



From (5.231, we can get the imarginary part (3 = ip as 

p 2 = (Pr + Lc + Pr Le ) 7 4 + Pr (1 + Le )(a 2 S 1 + ir 2 S ) + Pr a 2 n c /-/ 2 
= 7.3 x 10 16 . 
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Thus, we derive that 

^2 



i?2 



7T- 



Co 



2tt 
P 2 

Le 7T 2 (27T+ 1) 



#3 = 

C 4 = 



Le tt 2 

7r(7T 2 - 2) 



-Si = -s-(iZ — R C2 )i 



h=h = - 

C5 = — 

#2 



7r 2 (i? C2 - Le fl) 



Then the parameter b 2 in (|4.2|) reads 
(5.41) 



62 = -r [Le 7r(27r + l)i? 2 + (tt 2 - 2)i? 2 



P u 

Then Theorem 14.31 is rewritten as 



- (2tt 2 - tt - 2 + Le 7r(7r + 2))RR C2 ] 
r [_(l_ Le )^(7r-l)i? 2 + (7r 2 -2)^ 2 

+ (tt - 2 - Le vr 2 - 2Le ir)Rr) c ] 
':[-tt(tt - l)R 2 + (tt 2 - 2)rf c + Rr] c }. 



Physical Conclusion 5.3. Let R > 2.35 x 10 , and hi be the -parameter as in 



(5.4-1). Then for the problem (5.21) with (2.8), Assertions (1) and (2) of Theorem 
\4-3\ hold true. 

6. THC Dynamics 

The above Physical Conclusions 5.1-5.3 provide the possible dynamical behaviors 
for the great ocean conveyer, depending on the saline Rayleigh number R. 

We note that the temperature and salinity differences Tq — Ti and So — Si between 
the oceanic bottom and upper surfaces are different from the observed data. In fact, 
the observed values should be as follows 

AT = T Q - 7\ + T(x-l,x 2 ,0) - T(xi,x 2 , 1), 
AS* = So — Si + S(xi,x 2 ,0) - S(xi,x 2 , 1), 

where T{x) and S(x) are the transition solutions as discribed by Assertion (3) of 
Theorem |4.1 1 Due to the heat resources coming from the earth's crust, the bottom 
temperature To of the ocean retains essentially a constant, and the bottom salinity 
Sq is also a constant which equals to the average value 

So = — r / S(x)dx, \V\ the volume of the ocean. 
Im Jv 

It is easy to see that on the upper surface, the temperature and salinity vary in 
different regions. In the equator, T is about 20 ~ 40°C and near the Poles, T is 
2 ~ — 1.8°C. Here, T\ takes an average on the upper surface. Since the elevated 
rate of evaporation in the tropical areas and the freezing of polar sea ice, which 
leave the salt behind in the remaining seawater, the densities of tropical and polar 
water are quite high. Hence, as an average, we determine that 

(6.1) 



To - T x < 0, So - Si < 0. 
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By Physical Conclusion 5.1 and (6.1 ), the saline Rayleigh number R < 0. There- 
fore the transition of (5.211 with (2.8 1 is from real eigenvalues. In addition, by 
R < 0, we see that the number in (5.40) satisfies that b\ > 0. Thus, by Physical 
Conclusion 5.2, the dynamical behavior of the THC is a Type-I transition to a pair 
of stable cquilibrum states provided 

a = R 



Le 1 R > er c , 



which, by (5.271, (5.28) and (5.371, is equivalent to 

(6.2) 0.86 x 10 21 (T - T x ) + 3.75 x lO 23 ^ - S„) > 2.33 x 10 23 . 

For the large scale ocean circulation, T Q — Ti < 100°C. Hence, (6.2) shows that 
if Si — So = O(l), then the great ocean conveyer is driven by the doubly-diffusive 
convection. 

In Assertion (1) of Physical Conclusion 5.2, the steady state solutions describing 
the great Conveyer are appoximatively expressed by 

±C(3 1 / 2 < 



(6.3) 



(6.4) 



(6.5) 



(cr)L\ sin — — cos 7ra3, 



T ± = T 



=F Cj3 1 / 2 (a) cos ^p- sin ttx^ , 

(Ti - T a )x 3 =F — s-, — w- cos — - sin ttx 3 , 
od + 7T- 2 Li 



S 1 * = S + (Si - S )x 3 =F 



sign(5 - S^Cp 1 / 2 



TTXi . 

cos — — sm7ra;3, 



Le (a 2 + 7T 2 ) 

where C > is a constant, (3(a) is the first real eigenvalue, and by the eigenvalue 
crossing properties in |22j . (3(a) can be expressed as 

(6.6) P(o) = k(a - a c ) + o(\a - a c \) 

for some constant k > 0. 

In oceanic dynamics, the Rayleigh number a = R — Le _1 i? is a main driving 
force for the great ocean conveyer, (6.3 1 and (6.6 1 show that the velocity v of the 
oceanic circulation is proportional to y/a — a c . 

The velocity field v + given by (6.3 1 is topologically equivalent to the structure 



as shown in Figure [BTTj which is consistant with the real oceanic flow structure. 

Assertion (3) in the Physical Conclusion 5.2 shows also that the theoretical 
results arc in agreement with the real thermohaline ocean circulation. 

At x 3 — 1/2, the temperature T + and salinity S + given by (6.4 1 and (6.5 1 are 



as follows 



1 l/2 



1 



Tx) 

at 



TTXi 

-^cos — . 
r 2 L\ 



S+ /2 (x 1 )=-(S + S 1 ) + 



C/3 1 / 2 TTXl 

Le (a 2 + 7r 2 ) C0S T7' 



-i/2( x i) an( i ^i/2( Xl ) are illustrated by Figures 
spectively. In particular, if x\ — stands for the North Atlantic 



6.2 



and 



6.3 
Iot 



re- 
the 



The distributions T, 

xT= Li 

North Pacific, and £3 = 1/2 for the deep basin of the ocean, then the profiles of 
T^ 2 (x\) and S^ 2 (x\) confirm with observations, in where the North Atlantic Deep 
Water is cold and salty water, and the North Atlantic Deep Water is warmer and 
fresher water. 
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Figure 6.1. The flow structure described by v + in (6.3 1 




Figure 6.2. The temperature profile in deep water of the ocean. 




Figure 6.3. The salinity profile in deep water of the ocean. 



Physical Conclusions 5.1 and 5.2 show that although the THC is mainly due to 
a continuous transition to multiple equilibria, analysis for other possible dynamical 
behaviors is still interesting. 

We can see that when R satisfies 



(6.7) 



Le 3 - Le 2 

rfj c < R < ; — Or 



1 -he 2 



1 - Le 



24 



MA AND WANG 



the number b\ in (5.40) is less than zero, i.e. b\ < 0, then due to Physical Conclu- 



sions 5.1 and 5.2, the transition is jumping, leading to a saddle-node oscillation. 

When R > j^a c = 2.35 x 10 19 , the problem under goes a dynamic transition 
to a periodic solution, as rj > r/ c . This transition can be either TypeT, or Type-II, 
described as follows: 

First, if R > 2.35 x 10 23 , the number b 2 in (5.411 is less than zero, and by 



Physical Conclusion 5.3, this system undergoes a Type-I transition to a stable 
periodic solution as 77 = R — R > rj c , its velocity field is written as 

1/2 



(6. 



u{x, t) 



2ttA(?7) 

\ N J 



1/2 



L\ sm pot sm — — cos 7ra;3, 

■ t 7TX i ■ 
sm p cos — — sm ttx^ , 



where \{rf) — Re/3i(7/), 77 = R — R, and po is as in (5.41 I 



Second, if 2.35 x 10 19 < R < 2.35 x 10 23 , the number b 2 > 0. Hence the transition 
is jump to a stable periodic solution at 77 = 7/ c , leading to an oscillation between a 
time-periodic solution and the trival equilibrium near 77c . 

Finally, some remar ks fo r the oceanic thermohaline circulation are in order. 

First, the condition (6.1 I provides a more realistic scenario, where R satisfies that 
R < rzf-s o" c - In this case the circulation is a continuous transition to multiple 
stable equilibria, with flow structure as illustrated by Figure |6.1| and with the 
temperature and salinity profiles for deep water of the ocean as shown in Figures 



( 6.2 1-( 6.3 ) . This is consistent with observations. 



Second, it i s kn own that the velocity takes n/h as its unit. Thus, the maximal 



number given by (5.401, and 61 = 

kLi 



1/2 

value of vi in ( 6 . 3 1 is that v max — CLiK.hr 1 ft 1 / 2 , where C = (^p&i) , &i is the 



100 

8 



From (5.261 and (5.36) we get 



P 1/2 (o) = 0.64 x l(T 7 /? 1/2 (<7)m/s. 



The theoretical velocity is very small. In addition, the ratio between the vertical 
and horizontal velocities is given by 

— = — = 1.56 x 10 . 
vi Lt 

As a contrast, the oceanic circulation takes about 1600 years for its a journey, i.e., 
the real velocity is also very small. 



Third, the condition (6.7 1 would lead to a saddle-node bifurcation to metastable 



(6 



states. However, these state have not been observed in oceanography. 

Fourth, if R > 2.35 x 10 23 , then the circulation is time-periodic, and we see from 
that the period is 

2tt h 2 . 
T= = 1.1 x 10 6 s. 

Po « 

6.5 x 10 8 . This period is about four months. It 



where po is as in (5.41 1, and po 
is not realistic. 

Fifth, when 2.35 x 10 19 < R < 2.35 x 10 23 , the oceanic system undergoes a time- 
periodic oscillation phenomenon. However, this behavior has not been observed in 
a realistic oceanic regime. 
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7. Concluding Remarks 

Two criteria are derived in this article. First, a nondimensional parameter K is 
introduced to distinguish the multiple steady state and oscillatory spatiotemporal 
patterns. These patterns play an important role in the understanding the mecha- 
nism of thermohaline circulation in different oceanic basins. Second, for both the 
multiple equilibria and periodic solutions transitions, both Type-I (continuous) and 
Type-II (jump) transitions can occur, depending respectively on the signs of two 
computable nondimensional parameters b\ and 62- 

A convection scale law is introduced, providing a method to introduce proper 
friction terms in the model in order to derive the correct circulation length scale. 
The analysis of the model with the proper friction terms shows that the THC ap- 
pears to be associated with the continuous transitions to stable multiple equilibria. 

The study provides some general principles and methods for the dynamic tran- 
sitions and stability associated with thermohaline circulations, and will be used in 
different flow regimes in forthcoming articles. 

Appendix A. Dynamic Transition Theory for Nonlinear Systems 

In this appendix we recall some basic elements of the dynamic transition theory 
developed by the authors [171 [20], which are used to carry out the dynamic transition 
analysis in this article. 

In sciences, nonlinear dissipative systems are generally governed by partial dif- 
ferential equations, which can be put in the perspective of a dynamical system, 
finite or infinite dimensional, as follows: 

du 

(A.l) — = L x u + G(u,X), u(0) = uo, 

at 

where u : [0,oo) — ► H is the unknown function, A e I 1 is the system parameter, 
and if is a Banach space. We also need a Hilbert space Hi such that the inclusion 
Hi C H is compact and dense. 

Linear theory and principle of exchange of stability. Linear theory for 



system (A.l I is closely related to the principle of exchange of stability (PES), leading 
to precise information on linear unstable modes. To be precise, let {/3j(A) £ C | j £ 
N} be the eigenvalues (counting multiplicity) of L\, and assume that 

r <0 ifA<A , 

(A.2) Re /Jj(A) < = if A = A , VI < i < m, 

[ >0 ifA>A , 

(A.3) Re /3j(\ ) < Vj > m + 1. 

Much of the linear theory on stability and transitions is on establishing the PES. 
There are a vast literature devoted to linear theory including, among many others, 
[TJ [5] for classical fluid dynamics, and [53] for geophysical fluid dynamics. Some 
formulas on the derivatives of the eigenvalues with respect to the control parameter 
are derived in [22], and are used to verify the PES in a much easier fashion. 

Center manifold reduction. In many nonlinear problems, we need to reduce 
the infinite (or higher) dimensional system to a finite (or lower) dimensional sys- 
tem. The most natural way for this purpose is to project the underlying system 
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to the space generated by the most unstable modes, fully preserving the dynamic 
transition properties. This is achieved by using the center manifold reduction. 

To be precise, assume that L\ : Hi — > H is a parameterized linear completely 
continuous field depending continuously on A e I 1 , which satisfies 

L\ = —A + B\ a sectorial operator, 
(A. 4) A : H\ — > H a linear homeomorphism, 

B\ : Hi — > H a linear compact operator. 

In this case, we can define the fractional order spaces H a for a € K 1 . Then we also 
assume that G(-, A) : H a — > H is C r (r > 1) bounded mapping for some < a < 1, 
depending continuously on A € I 1 , and 

(A.5) G(u,X) = o(\\u\\ Ha ) VAeM 1 . 



Hereafter we always assume the conditions (A. 4 1 and (A.5 I, which represent that 
the system ( A.l I has a dissipative structure. Then the type of transitions for ( A.l ) 
at (0, Ao) is dictated by its reduction equation near A = A on the center manifold 
corresponding to the first m eigenvalues as given in (A. 2 1 and (A.3): 

dx 
~dt 

where g(x, A) = {gi{x, A), • • • ,g m (x, A)), and 



(A.6) 



J\x + g(x, A) 



for x £ 



(A.7) 



gj(x,X) =< G(Va; l e J + $(x,A),A),e* > VI < j < m. 



i=l 



Here ej and e* (1 < j < m) are the eigenvectors of L\ and L* x respectively cor- 
responding to the eigenvalues (3j(X) as in (A.2| and (A.3), J\ is the mx m order 
Jordan matrix corresponding to the first m eigenvalues, and $>{x, A) is the center 
manifold function of (A.l) near Aq. In addition, let 



H = Ei®E 2 , 
Ei = spanje^ 
£a = L\\e 2 - 



1 < i < m}, E 2 = E 



1 ! 



The center manifold function <!> is implicitly defined, and is oftentimes hard to 
compute. A systematic approach is developed in [15j[17l[20] to derive approxima- 
tions of 4>, which provide complete information on the dynamic transition of (A.6 1, 
consequently the original system (A.l I. Suppose the nonlinear operator G to be of 
the form 



(A.8) 



G(u, A) = G k (u, A) + o(||u|| fe ), as u -> in 



for some integer k > 2, where Gk is a k- multilinear operator 
G fe (it, A) = G k (u, ■ ■ ■ ,u, X) : Hi x ■ ■ ■ x Hi — 



H. 



Theorem A.l. [15] Under the conditions (A. 2), (A.3) and (A.8), the center man- 
ifold function <5>(x, A) can be expressed as 



(A.9) 



Pe P 2 G k (e rJ *x,X)dT + o(\\x\\*), 



where x = x^e^ € E\. In particular, we have the following assertions: 
i=i 
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(1) If J\ is diagonal near A = Ao, then (A.9\ ) can be written as 
(A.10) - £ x <S>(x, A) = P 2 G k (x, A) + o(\\x\\ k ) + 0(|/3||M| fc ), 

where /3(A) = (/3i(A), • • • , /3 m (A)) are the eigenvalues of J\. 

(2) Ietm = 2 and (3 2 {\) = (3 2 {\) = a(X) + ip(X) with p(A ) ^ 0. IfG k (u,X) = 
G 2 (u,X) is bilinear, then the center manifold function <&(x, A) can be ex- 
pressed as 

((-£ x ) 2 + 4p 2 (\))(-£ x )<S>(x,\) 

= ((-C x ) 2 + 4 P 2 (X))P 2 G 2 (x, A) - 2p 2 (\)P 2 G 2 (x, A) 

+ 2p 2 P 2 G 2 (xie 2 - x 2 ei) + p(-£\)[G 2 (xiei + x 2 e 2 ,x 2 ei - xie 2 ) 

+ G 2 (x 2 ei - xie 2 ,xie 1 + x 2 e 2 ) + o(||x|| 2 ) + 0(|a||ja;|| 2 ). 



(A.ll) 



Classification of dynamic phase transitions. A starting point of the dynamic 
transition theory is the introduction of a dynamic classification scheme of dynamic 
transitions, with which phase transitions, both equilibrium and non-equilibrium, are 
classified into three types: Type-I, Type-II and Type-Ill. Mathematically, Type- 
I, II and III transitions are also respectively called continuous, jump and mixed 
transitions. As we know, for equilibrium phase transitions, the usual classification 
scheme of phase transitions is based on the classical Ehrenfest classification scheme 
such that phase transitions are labeled by the lowest derivative of the free energy 
that is discontinuous at the transition. 

Here we give a brief description about this classification. A state of the system 
(A.l I at A is usually referred to a compact invariant set A state of ( A.l I is 
stable if E A is an attractor; otherwise T,\ is called unstable. 

The system (A.ll undergoes a dynamic phase transition from a state T,\ at 
A = Ao if is stable on A < Ao and is unstable on A > Ao. The critical parameter 
Ao is called a critical point. In other words, the phase transition corresponds to an 
exchange of stable states. 



Assume that we have the linear theory at our disposal, i.e., the conditions (A. 2) 



and (A. 3 1 hold true. Then we can show that the system (A.l I undergoes a dynamic 
transition from (u, A) = (0, Ao), and there is a neighborhood U C X of u = 
such that the transition is one of the three types, Type-I, II, and III as shown 
in Figures A.1|A.3| we refer interested readers to [20] : Type-I (continuous) 



Figure A.l. Schematic of Type-I transition: The transition 
states are represented by a local attractor £>, which attracts a 
neighborhood of the basic solution. 



transition is essentially determined by the attractor bifurcation theorem proved in 
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Figure A. 2. Schematic of Type-II transition: The transition 
states are represented by some local attractors which are away 
from the basic state at the critical A . 



Figure A. 3. Schematic of Type-Ill transition: The transition 
states are represented by two local attractors, with one as in a 
Type-II transition, and the other as in a Type-I transition. 



[151 [T7] . The attractor bifurcation theorem amounts to saying that when the PES 
holds true and the basic state is asymptotically stable at the critical parameter 
value Ao, the system undergoes a Type-I dynamic transition, which is described by 
the bifurcated attractor. The study of the attractor bifurcation theory was initiated 
a few years ago by the authors, and has been applied to many problems in sciences. 
The key assumption here is the asymptotic stability of the basic solution at the 
critical parameter value Ao- Two methods have been used to verify this condition 
in applications. One is a general alternative principle used in the study of Bernard 
convection |16l 117] . The other method is to use the center manifold reduction. 

When the asymptotic stability of the basic state at the critical parameter is no 
longer valid, the system undergoes either Type-II or III transitions, depending on 
the nonlinear terms. Hereafter we list a few theorems, which are used directly in 
proving the main results in this article. 



Transitions from simple eigenvalues. We consider the transition of ( A.l I from 



a simple critical eigenvalue. Let the eigenvalues (3j{\) of L\ satisfy (A. 2) and (A. 3) 
with m = 1. Then the first eigenvalue /3i(A) must be a real eigenvalue. Let ei(A) 
and e*(A) be the eigenvectors of L\ and L* x respectively corresponding to /3i(A) 
with 

£A o ei=0, L* x el = 0, <ei,ej>=l. 



Let Q(x, A) be the center manifold function of (A.l I near A = Ao- We assume that 
(A. 12) < G(xei + ®(x, Ao), Ao), ej >= ax k + o(\x\ k ), 

where k > 2 an integer and a/Oa real number. 
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(a) 



(b) 



Figure A. 4. Topological structure of the jump transition of ( A.l I 
when A:=odd and a > 0: (a) A < A ; (b) A > A . Here the 
horizontal line represents the center manifold. 



t 

(b) 



Figure A. 5. Topological structure of the continuous transition of 



(A.l I when k=odd and a < 0: (a) A < A ; (b) A > A . 



Theorem A. 2. Assume (A. 2) and \A.3) with m = 1, and {A. 12). If k=odd and 
a =/= in (A. 12) then the following assertions hold true: 



(1) If a > 0, then (A.l) has a jump transition from (0, Ao), and bifurcates on 
A < Ao to exactly two saddle points Vi and v£ with the Morse index one, 

H3 



as shown in Figure 



(2) lfa<0, then (A.l) has a continuous transition from (0,Aq), which is an 



attractor bifurcation as shown in Figure A. 5 
(3) The bifurcated singular points Vi and v 2 in the above cases can be expressed 
in the following form 

vl 2 = ±\p 1 (X)/a\ 1 '^e 1 (\) + oflftlVC*-!)). 

When /eleven and a ^ 0, one can prove that there will be a mixed transition, 
we refer the interested readers to [2U] for more details. 



Complex simple eigenvalues. We now study the transition from a pair of com- 
plex eigenvalues. Assume that the eigenvalues of L\ satisfy 



(A.13) 



(A.14) 



Reft (A) = Reft (A) 



<0 ifA<A , 
= ifA = A , 
> if A > A , 
I Imft(Ao) = -Imft(Ao) ± 0, 
Reft(A )<0 Vj>3. 
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It is known that with the conditions (A. 13 1 and (A. 14), (A.l I undergoes a Hopf 
bifurcation from (0,Ao). The following theorem amounts to saying that the tran- 
sition of (A.l I from (0, Ao) has only two types: continuous and jump transitions, 
which can be determined by the sign of a number b defined by ( A.18[ ). 

For this purpose, let ei(A),e2(A) and e*(A) and e 2 (X) be the eigenvectors of L\ 
and L* x respectively corresponding to the complex eigenvalues /^(A) = /3i(A) = 
a(A) + iS(\), where a satisfies (A. 13 1, a = cr(X ) ^ 0, and 

f L A e 1 (A)=a(A)e 1 (A) + <7(A)e 2 (A), 
\ L A e 2 (A) = -a(A) ei (A) + a(A)e 2 (A), 

( L* x el(X) = a(A)e* - a(X)e* 2 , 
\L* x e* 2 {\) = a(X)el + a(\)e* 2 . 

By the spectral theorem in |17j.we can take 

<e 4 (A),e*(A) >= VL<i,j<2. 



(A.15) 
(A.16) 



Let <fr(x, A) be the center manifold function of (A.l I near A = Ao, x = xiei+x 2 e 2 , 
and ej = ej(Ao)(« = 1, 2). Assume that for i = 1,2, 

(A.17) <G(x + $(i,Ao),A ) ) cJ>= £ o*,^ + o(\x\ 3 ). 

2<P+9<3 



For (A.17 1 we introduce a number, which is the bifurcation number: 
, 3tt 



(A.18) 



Here 



4 -(«30 + «03) + 4 Ol2 + a 2l) + ^( a 02 a ; 



1 2 \ 
a 20 a 2(V 



4ff 



(«ll a 2C 



_ 2 2 \ 
a ll a 02J- 



a(A ). 



Theorem A. 3. Lei £/ie conditions (A. 13) and (A. 14) hold true. 

(1) If b < 0, then the transition of {A.l) is continuous, and the bifurcated 
periodic orbit is an attractor. 

(2) lfb>0, then the transition is jump, and {A.l) bifurcates on A < Ao to a 
unique unstable periodic orbit. 



Singular separation. We now study an important problem associated with the 



discontinuous transition of (A.l I, which we call the singular separation. 



Definition A.l. 



(1) An invariant set £ of (A.l) is called a singular element 



if £ is either a singular point or a periodic orbit. 



(2) Let £i C X be a singular element of (A.l) and U C X a neighborhood of 



We say that {A.l) has a singular separation of £ at A = Ai if 

(A.l) has no singular elements in U as X < Ai (or X > Xi), and 



£i- 

(a) 

generates a singular element £i C U at X = X\, and 
(b) there are branches of singular elements £>, which are separated from 
£i for A > Ai (or X < Xi), i.e., 



lim max dist(x.Y,\) = 0. 
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Figure A. 6. Saddle-node bifurcation. 




Figure A. 7. Singular separation 
of periodic orbits. 



A special case of singular separation is the saddle-node bifurcation. Intuitively, a 



saddle-node bifurcation is schematically shown as in Figure A. 6 where the singular 
points in Ti(A) are saddle points and in r 2 (A) are nodes, and the singular separation 
of periodic orbits is as in shown Figure |A.7| 

The following theorem gives a general principle for singular separation. 



Theorem A. 4. Assume (A. 2) and (A.Sy . Then the following assertions hold true: 
(1) // (A.l) bifurcates from (u, A) = (0, Ao) a branch E> of singular elements 



on A < Ao which is bounded in X x (— oo,Ao) then \A.l ) has a singular 
separation of singular elements at some (Eo,Ai) C X x (— oo,Ao). 
(2) If the bifurcated branch T.\ consists of singular points which has index —1, 
i.e., 

ind(-(L x + G),u\) = -1 Vu A € E x , A < A , 

then the singular separation is a saddle-node bifurcation from some (u\, X\) G 
X x (—oo, Ao). 

Consider the case where /?i(A) = $2 (A) are a pair of complex eigenvalues of L\, 
and assume that 

< G(u, X),u >h= 

(A.19) 



< L\*u,u >#< — a \\u\ 



Vm e Hi, 
for some A* < Ao, 



(A.20) 



< G(u,\),u > H < -ci||u||^ + c 2 ||m||^, 

< L x u,u > H < ~a\\u\\ 2 H +cx\\u\\ 2 H , 



where a, C\, c% > are constants, p > 2, c\* < — c 2 for some A* < Aq, and Aq is as 



in the PES (A.13I and (A.14I 



Theorem A. 5. Under the PES {A. 13) and (A. 14-), the equation (A.l) has the Hopf 



bifurcation at (0, Ao). If the branch T.\ of bifurcated periodic orbits is on A < Ao ; 
then we have the following assertions: 



(1) // (A.19) holds, and for any A, (A.l) possesses a global attractor in H. 
then ( p4.jp has a singular separation of periodic orbits at some (Si,Ai) C 
Hx (A*,A ). 
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(2) If \ A. 20) holds true, then (A.ly has a singular separation of periodic orbits 
at some (Si, Ai) C H x (A*, Ao). 

(3) The branch Y.\ of bifurcated periodic orbits converges to Si as A — > Ai. 

Appendix B. Proof of the main theorems 



B.l. Proof of Theorem |3.1[ The proof of Theorem |3.1| is achieved using the 
following two lemmas. 

The the first lemma ensures that the R-Rayleigh number a defined by (3.151 is 
a reasonable parameter describing the critical-crossing for the real eigenvalues of 

Lemma B.l. Assume that 



-Le 2 (l + Pr) 7 6 ifcil 
(1 - Le )Pr a\ ki ' 



and for the R-Rayleigh number a near a c , all real eigenvalues of (3.1) are given by 

(B.l) 01>02>~->Pm> Pm+1 >■■ ■ 

Then (3\ — Ph^ii an d 

!< if a < a c , 

= if a = a c , 

> if a > cr c , 



(B.2) 



where a c is as in (1.2). 



Proof. Let a = a jlkl ,j = j jlkll ,Ro - Lc 1 R = 7 6 /a 2 , where (ii,fci) as in |l.2 l. 
We shall show that 



(B.3) 



(Pr + Lc + Pr Le ) 7 4 - Pr a 2 -/- 2 (R - R ) > 0. 



Assume that (B.3) is not true, we consider the case 

(Pr + Le + Pr Lc ) 7 4 - Pr a 2 7" 2 (i? - ^o) < 0. 



(B.4) 

Note that the solution (3(a) — flj 1 k 1 i( (J ) of (3.141 with (j,k,l) = (j'i,fci,l) is 
tinuous on cr, and /3(tr c ) = 0. Hence 

(B.5) (3(a) -> as a -> a c = R - Lc^Rq. 



Thus, near a — a c the equation (|3.14| can be written as 
(B.6) 
where 



P(«) = ^4 + o((3(a)), 



6i(<r) 



2 / 7 



(B.7) 



bo(a) = Pr Le a 



h(a) = (Pr + Lc + Pr Le ) 7 4 - Pr a 2 j 2 (R - R). 



It follows from (|B.4[)-(|B.6| that 

> if a < cr c 

(B.8) (3(a) { =0 if a = a c 

> if a > o> 
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for a = R — Le R near a r — Rn — Le 1 Rn 



(B.9) 



We write the equation (3.14) in the following form 

(3 3 + b 2 p 2 + b 1 (a)/3 + b (a) = 



where bo(a), &i(cr) are as in (B.7l, and 

b 2 = (Pr +Lc + 1) 7 2 > 0. 
Meanwhile, it is easy to see that 

b (a) — > +oo, bi(a) — > +oo, as u — 



It implies that when tr is sufficiently small, the real solutions of (B.9) must be 
negative. Thus, by (B.8) there exists a number do < o~ c such that the solution /3(a) 
of (B.9 1 vanishes at a — oq which is a contradiction to that &o(fo) 7^ 0- Thus, we 
derive b\(a c ) > 0. By the assumption in the lemma, bi(a c ) ^ 0. Hence, (B.3l holds 
true, i.e., &i(cr c ) > 0. By (B.6l, we can obtain (B.2|. 

In the following, we shall prove that ffi fci 1 — (3\ as in (B.l I. We only need to 
consider the real eigenvalues j3 satisfying (3.141. Let j3 m — Pjki( ~c) be a solution of 
(3.141 at er = er c . We consider the coefficients of (3.141 at a = <j c . Thanks to (1.2 1, 

(Pr + Le + 1) 7 ? H > 0, 



Pr Le - Pr a 2 fe (Le R - R) = Pr a%Le 



~6 



> 0. 



Thanks to (B.3l, we have 
(Pr 



Le + Pr Le h 6 



Pr a 2 



> Rq — Rq. 



Thus, we obtain 



(Pr + Le + Pr Le ) 7 * w - Pr a%^{Ro 

„6 



i2o) 



Pr +Le 



Pr 



Pr Lc 7^ 
7^2 \ K o 



Bo) 



> a %7 j kii Pl +Le +PrLe) 



> (by (|OJ) 



Hence, the coefficients of (|3.14|) at a = a c are nonnegative, and strictly positive 



provided that l%i/ot 2 kl ^ a c . It follows that 



/3 m = PjkMc) < V<r c ? 7 ? fcI /a 2 



j'fcl- 



Thus, we derive that (3j 1 k 1 i = /3i for <r near er c . The proof is complete. 



□ 



The following lemma shows that the C-Raylcigh number rj defined by (3.17) 
characterizes the critical-crossing at rj c for the complex eigenvalues of (3.1 ). 

Lemma B.2. Let (ji,ki) satisfy (3.19), and the condition 13.18 ) hold true for 
(j, k, I) = (ji, k\, 1). Then the pair of complex eigenvalues Pj^iiv) = Pj^xiv) are 
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critical-crossing at r\ = r/ c 



<0 ifr}<T) c , 

Ref3j lkll ( V ) \ =0 ifv = Vc 
>0 ifi]>r] c . 



Proof. Near r\ — r\ c the solution (3^ ikil {rj) of (B.9l takes the form 



(B.10) 



\(r)) -> 0, p{rf) -> po as-q^ri c 



Inserting (B.lOl into (B.9l we get 

(-3p 2 + 6i)A + b - b 2P 2 + o(A) = 0, 



(B.ll) 



- p 3 + pbx + 2pb 2 X + o(A) = 0. 



Since po ^ 0, we derive from (B.ll I that 

ba - hp 2 _b a - b x b 2 - 2b\X 



\{ri)+o{\) = 



3p 2 - h 



which yields 
(B.12) 
Note that 



2fei + 66 2 A 
b - b x b 2 



>(A), 



^ = 6i(r? c ) >0, A(r? c )=0. 



We derive from (|B.12|) that 
(B.13) 



< if fo < bih, 



Re$ lkl M = Hv) \ =0 ifb = b 1 b 2 , 
> if b > bib 2 , 

for 7/ near rj c . It is easy to check that 

< b x b 2 if 77 < ?/ c , 

fro \ = & i & 2 if r) = ?7c, 

> 6162 if V > rj c . 



Thus, the lemma follows from (B.13 1 



Proof of Theorem 3.1 Let a = R — Lc 1 R be at the critical state 
(B.14) R-Le~ 1 R = a c = j 6 /a 2 . 

Assume that 



(B.15) 



_ Pr +Lc - (Pr +Le )(1+Le ) 7 6 

v = R - , 1 R > Vc = — 



Pr +1 



Then, we deduce from (B.14I and (B.15) that 

Le 2 (Pr + l) 7 6 
(1 - Lc ) Pr «2 
Le 2 (Pr + l) 7 6 
(1 - Le ) Pr 



7? 



> 



< 



Pr 

if Le < 1, 
if Lc > 1, 



which implies that K < 0. 



□ 
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In addition, let r\ = t] c , namely 
Pr 



(B.16) 



R- 



Pr +1 



Le - (Pr +Le )(1 + Le ) 7 6 
K = 



Pi- 



Then noticing that K < 0, we can infer from (B.16 1 that 

(B.17) p 2 = (Pr + Lc + Pr Lc ) 7 4 - Pr a 2 ~f- 2 (R - R) > 0, 



Thus, by (B.17 1 and Lemma B.2 the conditions (B.14) and ( B. 15 1 imply that r\ c is 
the first critical Rayleigh number, and (3.22) and (3.231 hold true. 
Likewise, let 



(B.16 



r\ = R- 



Pr 



Pr 



Le - (Pr +Lc )(1 + Le ) 7 6 

R < 



1 



Pi- 



Then, it is clear that (B.14 1 and (B.18) imply that cr c is the first critical Rayleigh 
number, and 

Le 2 (Pr + 1) 7 6 



R 



(1 — Lc ) Pr a 1 
Le 2 (Pr + l) 7 6 
(1 — Le ) Pr a 2 



which implies > 0. 

By Lemma B.l the conclusions (3.201 and (3.211 are valid. Thus, the theorem 
is proved. □ 

B.2. Proof of Theorems |4.1j44.3l The proof of these two theorems is based on 
the dynamical transition theory briefly presented in the Appendix. The central 
gravity of the proof is to carry out the detailed calculation of the center manifold 
reduction of the original infinite dimensional dynamical system to a finite dimen- 
sional dynamical systems. 

Proof of Theorems \4-l\ and \4--S\ Let J-y = (ji, fci, 1), ipj 1 — ipj^i- The reduced 
equation of (2.6 (-(2.8 1 in H reads 

(B.19) ^=pj i{a)x + 



1 



where ip € H is written as 
(B.20) V = xi/jj, + 

$ is the center manifold function, and 
G(^i, $0 = -P((m • V)ua, (mi • V)T 2 , (ui • V)S 2 ) 



— (GW.V),^), 



(G(Vi,^),^s) 



0T 2 
uu-k—T 3 



dS 2q 

Miq ^3 

da;,- 



for ^ = (u i ,T <1 Si)e.ff(* = l,2,3). 

By Theorem |A.l | the center manifold function <i> satisfies that 



—Lx$ = x 2 P 2 G(ip Jl , V'jJ + high order terms. 
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Hence it is routine to calculate that 



(B.21) 

where 
and 



$ =x' 



^_ (Glib r. . M 



V2J + -8 1 (ib 1 ib 1 * ) ^ 002 

^0021^002, ^002) 



~ A)02 00002, ^002) 



V><1 



•o(2), 



02 j = 4Pr a 2 , -/^ = 4^ 2 , -/3 2 02 = 4Le tt 



"02J == i^2j — (fci<^2 sin 2j 1 ai7rx 1 cos 2fc 1 a27ra;2, 

— jiai cos 2jia\iTXi sin 2fcia27ra2, 0, 0, 0), 
0002 = 0002 = (0,0,0, sin 2^x3,0), 
^002=^002 -(0,0, 0,0, sin 2^3). 
Inserting (B.20I into ( B.19| , we have 

(B.22) 
where 



dx 



8(°c) = E 
1 



(U P\(&)x -&{&):,? + 0$), 

(g^ ^ ) [(G(^, 7l ,r Jl ),iPj) + ( g ^j d^h ^ 

{Gjibj, m, ^jKg^j, M 

-l32.j(lp2J,lb*j) 



(G(ipj 1 ,4>Ji)i ^002 ) (G(ib Jl , ib* Ji ) , ^002) 

— 1^002 (V'002, V'002) 
(Giibj, ,ib j, ) , V 2 o2 ) (G( ^ ) , ^02) 



/3o02 (V , 002 ) 0002 ) 

-- r {kia 2 u Jl - Jia 1 vj l )(kia 2 u* 1i - Jxotxv*^) T^T^ S Jl S* Ji - 



Here 



Pr a 4 



r 4 Le 



jiai7r 



a- 



Sj 



sign(5 - Si) 
Le 7 2 + /3} 1 ' 



n = 



S 



PtR 

Pri? sign(5 - Si) 



■/i 



Le 7 2 + /3 1 



By Pj 1 (a c ) = 0, direct calculation yields that at a = a c 



5(cr c ) 



Pr(Lc 3 i? - R) 



Le J i? - R 



87T 2 7 4 Le (tf>j 1 , ib*j 1 ) 8 ^2 Le 3 + ^ _ _S 



Putting R — Le ii + a c into <5(cr c ) , we derive that 



5 = S(a c ) 



1 Le 3 cr c - (1 - Lc 2 )^ 
87r 2 Le 3 (l+i)a c -^i?' 
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Hence using K > 0, the sign of S is the same as the sign of b\ defined by ( 1.3 1, and 
then the theorem follows from Theorems A. 2 and A. 4 and from (B.22I. □ 

Proof of Theorem \4~3\ L et J = (ji,Q, 1) and (i(rj) = /3j(rj) = A(r?) + ip{rj) be the 
first eigenvalue of (B.lOl near r\ — r\ c with 

(B.23) A(r? c ) = 0, p(r) c ) = p > 0. 

By (B.23 1, the eigenvectors tpj = ipj + iipj corresponding to f3(j]) are given by 



1 



■ sin <f> cos 77x3 , , cos </> sin 77x3 , 



ReA2(f3) cos <f> sin 77x3 , Re A3 {(3) cos 4> sin 77x3) , 
(0, 0, 0, ImA2(/3) cos <fi sin 77x3 , Im^ cos <f> sin 773:3), 



where <ft 
given by 



iij} 2 j* are 



jiai77Xi. By (3.131 the conjugate eigenvectors ipj 

( sin </> cos 77x3, 0, cos 0! sin 77x3, 

Pr RReA2((3) cos </> sin 77x3, — Pr RRcA 3 ((3) cos </> sin 77x3), 
= (0,0,0, Pr i?ImA 2 (/3) cos sin 77x3, —Pr i?Iim4 3 (/3) cos <j> sin 77x3) 



and 
(B.24) 



a 2 + 7r 2 + /3' A3(/3) Le (a 2 + 77 2 ) + /5' 
The conjugate eigenvectors <&j* and <I> 2 *, satisfying 

= (^,$ 2 /) ^ 0, (^,$ 2 /) = (^,^*) - 0, 

are given by 



(B.25) 
where 

(B.26) 



iff, 



C = 



The reduced equations of (|2.6|)-(|2.8|) read 
(B.27) 



(ft 



= Ax + + 



^ (G(^) ,<!>?), 



where ^ € is as 
(B.28) 



and $ is the center manifold function. 
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Note that for any gradient field Vip, the Leray projection P(V<p) = 0. Therefore, 
by ( |R^ -( |R26) we have 

G(ip},ipj) = ^(0,0,0,ReA 2 sin27ra;3,ReA3sin27rx3), 
G(?Aj,V'j) = ^ (0, 0, JmA 2 sin 2^2:3 JmA 3 sin 27r£ 3 ), 



TrPr 



(0, 0, 0, R{B{ReA 2 - B 2 lmA 2 ) sin 2ttx 3 



R(B 1 ReA 3 - B 2 ImA 3 ) sin 2ttx 3 ) 



TrPr 
2S7 



(0, 0, 0, ~R(B 2 ReA 2 + Bilim4 2 ) sin 2nx 3 , 



R(B 2 ReA 3 + BihaA 3 ) sin 2712:3) 
G(^,V) = 0, V^etf. 



Hence, by (A. 11 1, the center manifold function $ is expressed as 
(B.29) 



$ = $ 1 +$ 2 + $3 + o(2), 



where dj — G(ipj,ipj), and 



/^002ll ' DII2 



2 112 I 



$ 2 = 



^ 2 [G 11 ,V4 2 ) + ^(Gl2,V'002)] i 



[(x 2 - y 2 )(G Ul 4 2 ) + 2xy(G 12 , ^02)] 



[(x 2 - 2/ 2 )(G n , V4 2 ) + 2zt/(G 12 ,^ 02 )], 



/3o 1 02(/5o 1 2 2 + V)IIV4 2 l| 2 

2/^002 

/3o 2 02((/3o 2 2) 2 + V)ll^0 2 02l| 2 

((^02) 2 +t 2 )li^2ll 2 [ ^ (Gll '^ ) + ^ - ^ 2)(Gl2 '^ )] 



002 



N(Gn,4) + (y 2 - z 2 )(G u ,</> 2 02 )]. 



((/3o 2 02) 2 + V)ll^o 2 o 2 ll 2 ' 
It is clear that 

(G($,$S*),^)=o(2) V*,j = l,2. 
Then, inserting ( B.28[ ) into ( |B.27 ), one gets 

Tt- Xx+py WW) + o(3) ' 



(B.30) 



a = + Ay - w>$a/) + o(3). 



From QB.29| ) and ( |B.30| it follows that 

Pr irxlaix 2 + a 2 xy + a 3 y 2 ] 

- „x + py + 

(B.31) 



dx 
~dl 
dy 
dt 



o(3), 



^j,^/) 2 + ^j^f) 2 

v Pr 7ra;[6ia: 2 + 6 2 xy + fe 3 y 2 ] 
= -px + Ay + ^ ,,i M2j _, ,1 ,2^2 + °( 3 )> 
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where 



01 = R{-R 2 B 1 + I 2 B 2 )D 1 + R(R 3 B 1 - I 3 B 2 )F 1 
a 3 = R{-R 2 B 1 + I 2 B 2 )D 3 + R(R 3 B 1 - I 3 B 2 )F 3 
b 2 = R{I 2 B l + R 2 B 2 )D 2 - R(I 3 Bi + R 3 B 2 )F 2 



where Bi,B 2 are as in (B.26|), and 
R 2 



Di 
D 2 
D 3 
Fx -- 
F 2 -- 
F 3 -- 



P 2 R 2 



(ml 2 



Then we obtain 
(B.32) 



8tt 167r(47r 4 + p 2 ) 4(4tt 4 + / 9 2 )' 
h p 2 h PttR 2 



in 8vr(47r 4 + p 2 ) 4(4tt 4 + p 2 ) ' 
h 



p 2 R 2 pirl 2 



167r(47r 4 + p 2 ) 8(4tt 4 + p 2 ) : 
R 3 p 2 R 3 



ph 



8?rLe 167rLe (47r 4 Le 2 + p 2 ) 8(47r 4 Le 2 + p 2 ) ' 

_h fy% pR 3 

87rLe 87rLe (47r 4 Le 2 + p 2 ) 4(4tt 4 Lc 2 + p 2 ) ' 

P 2 R 3 Ph 

167rLe (47r 4 Le 2 + p 2 ) 8(47r 4 Le 2 + p 2 ) ' 



b = 3ai + a 3 + b 2 

= |[i?(i? 2 S 1 - I 2 B 2 )C 3 + RihBx + R 2 B 2 )C A 
+R(R 3 B 1 - I 3 B 2 )C 5 + R(I 3 B X + R 3 B 2 )C e ] 



where b 2 is the number defined by (4.2 1. 

By (B.31 )-(B.32 1, this theorem follows from Theorems A. 3 and A. 5 The proof 
is complete. □ 
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